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Results from kinematic solar dynamo models employing a-effect and turbulent pumping from local convection calcula- 
tions are presented. We estimate the magnitude of these effects to be around 2-3 m s~ , having scaled the local quantities 
with the convective velocity at the bottom of the convection zone from a solar mixing-length model. Rotation profile of 
the Sun as obtained from helioseismology is applied in the models; we also investigate the effects of the observed surface 
shear layer on the dynamo solutions. With these choices of the small- and large-scale velocity fields, we obtain estimate of 
the ratio of the the two induction effects, C a /Cn w 10~ 3 , which we keep fixed in all models. We also include a one-cell 
meridional circulation pattern having a magnitude of 10-20 ms~ near the surface and l-2ms _1 at the bottom of the 
convection zone. The model essentially represents a distributed turbulent dynamo, as the a-effect is nonzero throughout 
the convection zone, although it concentrates near the bottom of the convection zone obtaining a maximum around 30° of 
latitude. Turbulent pumping of the mean fields is predominantly down- and equatorward. The anisotropies in the turbulent 
diffusivity are neglected apart from the fact that the diffusivity is significantly reduced in the overshoot region. We find 
that, when all these effects are included in the model, it is possible to correctly reproduce many features of the solar activity 
cycle, namely the correct equatorward migration at low latitudes and the polar branch at high latitudes, and the observed 
negative sign of B r B^. Although the activity clearly shifts towards the equator in comparison to previous models due 
to the combined action of the a-effect peaking at midlatitudes, meridional circulation and latitudinal pumping, most of 
the activity still occurs at too high latitudes (between 5° . . . 60°). Other problems include the relatively narrow parameter 
space within which the preferred solution is dipolar (AO), and the somewhat too short cycle lengths of the solar-type solu- 
tions. The role of the surface shear layer is found to be important only in the case where the a-effect has an appreciable 
magnitude near the surface. 
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1 Introduction 



The early solar dynamo models, initiated by Parker (1955), 
relied on a scalar a-effect and a negative gradient of the an- 
gular velocity, f2, in the convection zone (e.g. Steenbeck & 
Krause 1969; Deinzer & Stix 1971; Kohler 1973). These 
simple models were enormously succesful in reproducing 
many of the observed features of the solar cycle, most 
prominently the equatorward propagating activity belts. 

During the recent years, however, helioseismology has 
revealed that the radial differential rotation shows a positive 
gradient near the equator, whilst a negative gradient occurs 
only at high latitudes. Moreover, the strongest gradients oc- 
cur near the poles and at the bottom of the convection zone 
(e.g. Schou et al. 1998; Thompson et al. 2003). The sign 
of the J7-gradient, if combined with positive (negative) a- 
effect in the Northern (Southern) hemisphere, leads to pole- 
ward migration of the activity belts near the equator and 
to equatorward migration in the polar regions according to 
mean-field dynamo theory (Parker 1955; Yoshimura 1975). 
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This is in contradiction with the observed migration of solar 
activity tracers, being equatorward at low latitudes, whilst 
a weaker polarward wave occurs at high latitudes. The ob- 
served strong differential rotation in the polar regions, to- 
gether with the commonly adopted description of the a- 
effect being proportional to cos 8 therefore peaking at the 
poles, leads to dynamo solutions with high-latitude activ- 
ity; virtually no sunspots, however, are observed above the 
latitude ±40°. 

A suitable meridional circulation pattern is one possi- 
ble remedy to the aforementioned problems. At the mo- 
ment it is known that near the surface of the Sun, i.e. the 
topmost 10-20 Mm, there is a poleward flow of the order 
of 10-20 ms" 1 (Zhao & Kosovichev 2004; Komm et al. 

2004) , but very little is known about the structure of the 
flow in deeper layers. Models of the solar differential rota- 
tion predict that there should be a return flow near the bot- 
tom of the convection zone with a magnitude of approxi- 
mately 1ms" 1 (e.g. Kitchatinov & Riidiger 1999; Rempel 

2005) . The meridional flow structure and a low value of the 
turbulent diffusivity, rjt, are of critical importance for the so- 
called flux transport dynamos (e.g. Choudhuri, Schussler & 
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Dikpati 1995; Dikpati & Charbonneau 1999). These mod- 
els are essentially ail-dynamos where the meridional flow 
works as a conveyor belt for the poloidal magnetic field, 
and the a-effect is only due to the decay of active regions 
near the solar surface, or due to instabilities in the magnetic 
layer below the convection zone (e.g. Schmitt, Schiissler & 
Ferriz-Mas 1996; Brandenburg & Schmitt 1998; Dikpati & 
Gilman 2001), whereas the shear in the tachocline is respon- 
sible for the generation of the toroidal field. However, in 
the present paper we do not exclude the possibility to have 
an a-effect due to cyclonic turbulent motions in the con- 
vection zone proper as is expected in the mean-field the- 
ory and suggested by the local convection calculations (e.g. 
Brandenburg et al. 1990; Ossendrijver, Stix & Brandenburg 
2001, Ossendrijver et al. 2002; Kiipyla et al. 2006, hereafter 
KKOS). In this case it is questionable whether the merid- 
ional flow alone is able to counter the missing activity and 
wrong migration of the activity belts near the equator. 

Another possibility to solve the problems outlined above 
arises in the form of dynamo coefficients. Most of the so- 
lar models have up to date relied on a scalar a-effect, 
or they have taken into account a limited amount of the 
anisotropies, such as vertical pumping of the mean magnetic 
field (e.g. Brandenburg, Moss & Tuominen 1992), or the 
contributions to the a-effect due to the steeply decreasing 
turbulence intensity towards the bottom of the convection 
zone (Riidiger & Brandenburg 1995). Recently, Ossendri- 
jver et al. (2002) computed all of the coefficients oy that 
appear in the expansion 

— - + • • ■ 

dxk 

describing the small-scale effects on the large scales from 
local numerical 3D models of stellar convection. Taking 
only the first term on the rhs into account, £ can be writ- 
ten as 



Si 



bijk- 



(1) 



£ = olB + 7 x B 



where a. and 7 are defined via 

' a ji) > 



1, 



(2) 

(3) 
(4) 



and consist of the symmetric and anti-symmetric parts of 
fly. The diagonal components of ay describe the genera- 
tion of a large scale magnetic field through the a-effect. The 
vector 7 describes turbulent pumping of the mean field with 
a velocity that is common for all field components. The off- 
diagonal components of ay contribute to the field direction 
dependent part of the pumping effect (see e.g. Ossendrijver 
et al. 2002; KKOS). 

The main findings for the a-effect in the study of 
Ossendrijver et al. (2002) were that whereas it is highly 
anisotropic, the component corresponding to is almost 
constant in magnitude as a function of latitude for Coriolis 
number of Co = 2 fir « 2.5, where r is an estimate of 



the convective turnover time. It was also established numer- 
ically that the turbulent pumping of the mean magnetic field 
depends upon the field component. The results also clearly 
indicate down- and equatorward pumping of the mean mag- 
netic fields that can alleviate the problems facing solar dy- 
namo models. Furthermore, recently KKOS extended the 
study of Ossendrijver et al. (2002) to the rapid rotation 
regime, i.e. Co ~ 10, corresponding to the deep layers of 
the solar convection zone. The main result of this study is 
that a^^, responsible for the generation of the poloical field 
from the toroidal one, no longer peaks at the poles but rather 
at around latitude 9 = 30°, which could be of further help 
for the solar dynamo models employing the helioseismic ro- 
tation profile. The vertical pumping of the magnetic field is 
downward near the poles, but can be upward near the equa- 
tor. The latitudinal pumping is always towards the equator, 
but it is concentrated in a rather narrow latitude range, i.e. 
|0| < 30° near the equator. The vertical pumping of the 
toroidal field was found to be very small or upward at lati- 
tudes |6| < 45°. 

In the present study we explore the implications of the 
local 3D results for the solar dynamo by means of axisym- 
metric mean-field dynamo models employing the observed 
internal rotation of the Sun. Starting from a simple a-profile 
(cos 9 in latitude and constant in radius within the convec- 
tion zone) we demonstrate the problems facing solar dy- 
namo models. Then we introduce an a-profile that captures 
the essentials of the magnitude and latitude dependence 
found in local convection calculations. Then we proceed by 
adding the pumping effects and finally a reasonable merid- 
ional flow constrained by the flow velocity observed on the 
solar surface. We also discuss the stability of quadrupolar 
and dipolar field configurations for the chosen £1- and a- 
pro files with and without turbulent pumping and meridional 
flow. Finally, we discuss briefly the phase dilemma of the 
radial and toroidal magnetic fields which has been widely 
used as a constraint for solar dynamo models. The remain- 
der of the paper is organised as follows: in Sect.|2]the mean- 
field model is described in detail, and in Sects. [3]and|4]we 
present the results and the conclusions of the study. 

2 The model 

We study axisymmetric kinematic mean-field models of the 
solar dynamo in spherical polar coordinates within a shell 
O.6i?0 < r < Rq. The inner region of the shell up to r c = 
O.7i?0 models an overshoot region below the convection 
zone. We solve the mean-field induction equation 

— = V x[(U + ~f)xB + aB-r lt V xB], (5) 

where U = ft x r + (7 mer is the velocity, r — re r the 
radius vector, 17 = f2(r, 6)k the prescribed rotation pro- 
file, k the unit vector along the rotation axis, and [/ mcr = 
(U r ,Ug,Q) the prescribed meridional flow, a and 7 are 
given by Eqs. l|3} and (0), respectively. r]t is the turbulent 
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diffusion which we treat as a scalar field neglecting its ten- 
sorial nature for the time being. However, we take into ac- 
count the decrease of rjt in the overshoot region via 



Vt(r) =rj c + 0.4577o 



1 + tanh 



di 



(6) 



where r\ c — O.I770, d\ = 0.015, and 770 the value in the 
convection zone. 

Under the assumption of axisymmetry, the magnetic 
field can be represented with two scalar fields A and B ac- 
cording to 

B = B T + Bp=Vx (Ae*) + Be . (7) 

Thus we obtain separate evolution equations for A and B, 
which we non-dimensionalize by choosing the following 
units: 

[r] =R Q = R, [t} = R 2 / m = T, [U} = r] /R, 

[r) t ] = no, [B] = B , [A] = RB , [Q] = Q , (8) 

where the subscript refers to a typical value of the quantity 
in question. The resulting equations are given by Eqs. JA1I 
and iA2\ in Appendix A. The models are controlled by the 
dimensionless parameters 

a R fl R 2 ~ u R 



C a — 



Cn = 



C„ = 



(9) 



Vo Vo Vo 

describing the magnitude of the a-effect, differential rota- 
tion and meridional flow with respect to diffusion. In what 
follows, cto and uq are the maximum values of the a-effect 
and the latitudinal component of the meridional flow. 

The nonlinear back reaction of the magnetic field on the 
turbulence and thus on the a-effect is modelled via a simple 
a-quenching formula 



a l3 [r,9,B{r, 9)} 



l + [B(r, 6)1 U 



(10) 



where B cq is the equipartition strength of the large-scale 
magnetic field. In practise the value of B cq determines the 
saturation strength of the magnetic field and we set simply 
B eci /Bo = 1. The actual equipartition value can be esti- 
mated from typical values of the local turbulent velocity and 
density using a mixing length model of the solar convection 
zone. The values u Ims » 10ms _1 and p = lOOkgm -3 
corresponding to the situation at the bottom of the convec- 
tion zone yield saturation strengths of the order of a few 
kG, and similar values are obtained throughout the convec- 
tion zone. The same quenching formula is also applied on 
the turbulent pumping effects. 

The equations JA1I and (IA2> are solved over the radius 
0.6R < r < R and colatitude < 9 < tt using a two- 
dimensional equidistant grid of N r x Ng gridpoints. The 
boundary condition at the poles reads A = B = 0, and in 
the radial direction pseudo-vacuum conditions 
d{rA) 



B 



dr 



0, 



(11) 



are used at the surface, and perfectly conducting condition 
at the bottom of the convection zone, i.e. 



.4 



d{rB) 
dr 



Otggr 



OA 
dr 







(12) 



The numerical resolution is varied between 41 x 81 and 
61 x 121 grid points. The constant timesteps used are At = 
2.5- 10 -5 for the low resolution case and At = 1.5-10 _5 for 
the higher resolution case. The code has been validated us- 
ing the "dynamo benchmark" test cases (Arlt et al. 2006); 
the pseudo-vacuum boundary conditions used here were 
found to result in slightly reduced critical dynamo numbers 
in the case of ail dynamos if compared to methods with 
real vacuum boundaries, but otherwise good agreement with 
other methods was found. Due to this we prefer to keep the 
pseudo-vacuum boundaries at the surface due to the sim- 
plicity of their implementation. The code employs second 
order accurate explicit spatial discretisation and first (Eu- 
ler) or second order (Adams-Bashforth-Moulton predictor- 
corrector; see e.g. Caunt & Korpi 2001) accurate time step- 
ping. 

The cf-tensor and 7-vector component profiles are 
adapted from the local convection model of KKOS (see be- 
low). The Coriolis numbers realised in the convection mod- 
els are here interpreted as depths in the convection zone (see 
Fig. 1 of Kapyla et al. 2005), so that Co = 1 corresponds to a 
depth of roughly 50 Mm below the surface, Co = 4 matches 
the middle, and Co = 10 roughly the bottom of the solar con- 
vection zone. We neglect the possibility to have an a-effect 
in the overshoot region in the present study. The latitudinal 
dependences are approximative fits to the local calculations, 
made at four different latitudes for the sets Co = 1 and Co = 
4 with 30° steps, and seven different latitudes for the set Co 
= 10 with 15° steps. The rotation profile used throughout 
the investigation is a digitisation of the helioseismological 
inversion (e.g. Schou et al. 1998), see Fig.[0 Since no re- 
liable inversions are available for the polar regions, we as- 
sume the profiles of Q at |6| = 75° and |0| = 90° to be 
proportional to the profile at 60° via 



ft(Y,90 o ) = Q{r,60°) 



0.05 



rp-r 
r -R 



O.K^)-l 
ro-R 



-1 2 



(13) 
(14) 



where r = 0.6, and 6 = 90° - 9 is the latitude. With this 
choice fl at the surface on the poles is consistent with the 
observed surface value. 

In order to characterise the resulting solutions of the 
models in terms of the large scale magnetic field structure, 
we define the parity as 



P = 



E (S) + E (A) ' 



(15) 



where E^ and E^ are the energies of the symmetric and 
antisymmetric parts of the field, respectively (see, e.g. Bran- 
denburg et al. 1989). For a symmetric quadrupolar solution 
(SO) the parity is +1, and for an antisymmetric dipolar so- 
lution (AO), P = -1. 
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3 Results 

3.1 Estimates of C a and Co for the Sun 

Using the profiles of the transport coefficients from theoret- 
ical considerations or local 3D convection calculations, the 
remaining task is to choose appropriate values for the three 
dimensionless control parameters C a , Cq, and C\j, defined 
via Eq. l|9}. The local convection calculations of KKOS 
yield max(|ay |)/u rms w 0.2. We scale this to physical 
units by using the mixing length estimate of the turbulent 
velocity at the bottom of the convection zone, i.e. uJ^ T « 
10ms _1 , and obtain a magnitude of 2ms" 1 for cto. Fur- 
thermore, using the mean angular velocity of the Sun, fin = 
2.6 x 10~ 6 s~\ the ratio C a /Cn = ao/(0 -R) is therefore 
approximately 10 -3 ; in the dynamo models describing the 
Sun, we keep this ratio fixed. Furthermore, the meridional 
flow at the surface, uo » 10 . . . 20 m s~ 4 , is known from ob- 
servations, from which we obtain another useful ratio for the 
solar models, i.e. C v /Co = u /(Q R) « 0.005 . . . 0.01. 

The value of 770 still remains a free parameter. Sim- 
ple mixing length estimates for the deep layers of the so- 
lar convection zone (e.g. Stix 2002) give r)t = ^uf^d « 
5 • 10 8 m 2 s _1 , where d — ckMLT-Hp, «mlt = 1-66, and 
H p « 5 ■ 10 7 m. Similar or slightly higher values are ob- 
tained for all depths in the solar convection zone (see e.g. 
Brandenburg & Tuominen 1988). We find that growing so- 
lutions are possible when 770 ~ 10 8 m 2 s^ 1 (see below). 



Keeping the ratio C a /Co fixed at 10 -3 we find that the 
critical dynamo number, D CIlt = Qj, rlt Co, for the exci- 
tation of the antisymmetric AO mode is roughly £) cnt « 
4.89 • 10 4 corresponding to C£ rit w 6.99 using the att- 
approximation (D crlt w 4.37- 10 4 corresponding to C™* w 
6.61 if the a-terms are retained in the equation of B). In the 
marginal case the turbulent diffusivity turns out to be rjo w 
2.0 • 10 8 m 2 s . The toroidal magnetic field of the marginal 
solution exhibits an oscillation period of rj cyc « 0.044r, 
which corresponds to roughly 3.2 years in physical units. 

Considering the resulting butterfly diagrams we com- 
pare models where the shear above r = 0.95i? is either 
turned off (Fig. |5Jl or retained (Fig. In the former case 
the butterfly diagram of the azimuthal field near the sur- 
face shows that the migration of the activity belts at lati- 
tudes |0| < 60° is poleward with an equatorward branch 
near the pole, whereas the latter case exhibits a strong equa- 
torward propagating activity belt from latitude ss 80° all 
the way down to the equator. The migration of the former 
type has frequently been reported in mean-field models with 
distributed a-effect combined with helioseismological rota- 
tion profile (e.g. Bonanno et al. 2002). Concerning the latter 
case, we note here that if there is a non-vanishing a-effect 
near the surface, the surface shear layer plays a very impor- 
tant role in shaping the butterfly diagram. However, if the 
a-effect is more concentrated in the deep layers of the con- 
vection zone the surface shear plays only a minor role (see 
Sect.|33}. 



3.2 Constant a with cos 8 latitude dependence 



3.3 a-effect from the local convection models 



Although the main emphasis of this paper is to investigate 
the role of anisotropies found in the local convection mod- 
els, we perform a set of calculations with a constant a- 
effect in the convection zone with a previously commonly 
used cos 8 latitude profile in order to demonstrate the cur- 
rent problems of solar mean-field dynamo models. In all of 
the calculations in this section, we put 



1 



a, 



a Si 



1 + tanh 



di 



(16) 



where r c ~ 0.7 and di = 0.015. The tanh-factor takes 
care of the vanishing a-effect in the overshoot region. We 
study this simple case first and determine the critical dy- 
namo number C" 1 * as a function of Co for the antisymmet- 
ric AO and symmetric SO modes. The results for the AO case 
are shown in Fig. [2] We note that the critical values for the 
SO mode are very close to those of the AO mode, but consis- 
tently slightly larger. The difference, however, is less than 
one per cent in all cases. For values of Co larger than 10 4 
the difference between the aft and a 2 Vl models is less than 
5 per cent. In what follows we consider the case C a = 15, 
Co = 1.5 • 10 4 as our standard model, in which case we 
adopt the afi-approximation. Moreover, we note that the 
critical dynamo number in the case of no surface shear is 
roughly 20 per cent larger than in the case where it is taken 
into account. 



Already the early analytical investigations of the turbu- 
lent electromotive force showed that the a-effect is highly 
anisotropic (Steenbeck, Krause & Radler 1966). These 
anisotropies have recently been studied by means of local 
convection calculations by Ossendrijver et al. (2002) and 
KKOS. In the present study we employ these numerical re- 
sults of the a-effect and turbulent pumping in solar dynamo 
models. As was discussed earlier, we set the ratio C Q /Co to 
10~ 3 in the case of the solar dynamo. Effectively the small 
value of the ratio C Q /Co means that the dominant (diago- 
nal) component of a is which appears in the equation 
of the poloidal field (see Eqs. lAll and lA2> . Adopting the aVL- 
approximation and setting a rr = aee — does not signifi- 
cantly change our results. These components can, however, 
be of vital importance for a 2 -dynamos and their oscillation 
properties such as those investigated by Rudiger, Elstner & 
Ossendrijver (2003). 

In order to make use of the local results of KKOS we 
use the Coriolis number to determine the depth of the lo- 
cal convection model in the convection zone. In an earlier 
study we computed the Coriolis number from a mixing- 
length model of the solar convection zone (Fig. 1 of Kapyla 
et al. 2005), and found that it varies between 10 -3 near the 
surface and 10 or larger in the deep layers. In KKOS, sets of 
convection calculations at different latitudes are made with 
approximate Coriolis numbers of 1, 4, and 10. Identifying 
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Fig. 5 Left panel: profiles of the a^-component from three depths r = (0.765, 0.866, 0.946)i? corresponding to Coriolis 
numbers Co = (1,4, 10), respectively, using the convective velocity at the bottom of the convection zone as normalisation. 
Right panel: the actual a^r, 9) profile used in the models. 



these values with depths in the solar convection zone, they 
correspond to r = 0.946-R, r = 0.866i?, and r = 0.765R, 
respectively. In KKOS it was found that for Co = 1 the com- 
ponent corresponding to is consistent with a cos 9 lati- 
tude profile, whereas for Co = 4 is basically constant as 
function of latitude from the poles to the latitude |0| = 30°. 
In the rapid rotation regime, i.e. Co = 10, a.$$ peaks around 
|0| = 30° and the magnitude at the poles decreases. We use 
the following approximations of the latitude dependence of 
at the different depths 



t» w (r = 0.946i?,(9) = 
a<H> (r = 0.866i?,6>) = 



—a.Q cos 9 , 

— an tanh 
2 



'(»-!) 



(17) 



(18) 



a H {r = O.765.R,0) = 3 a a ( sin 2 9 cos 9 + - cos 6>) (19) 

see the left panel of Fig.[5](see also Figs. 3 to 7 of KKOS). 
We use cubic spline interpolation in order to obtain a value 
of for all grid points, see the right panel of Fig. |5] In 
Eqs. (I17> to dl9l an estimate of the convective velocity, 
u « 10ms _1 , at the bottom of the convection zone has 
been used to scale the profiles at all depths resulting an 
profile that peaks near the bottom of the convection zone at 
around latitude 30°. 

Considering models with the a-effect depicted in Fig. [5] 
we find that the critical dynamo number is significantly 
larger than in the case of constant a in the convection zone, 
see Fig. [6] Again the critical value for the SO mode is very 
sligthly larger than that of the AO mode. For the solar case of 
C a /C n = 10~ 3 we find £> crit w 8.2 • 10 4 , or C£ rit = 9.04, 
instead of D crit w 4.89 • 10 4 and C£ rit w 6.99 in the case 
of constant in the whole convection zone. This is of 
course due to the fact that the integral of over the full 
shell is significantly less than in the case of radially constant 
a.$$. Furthermore, we find that initially purely antisymmet- 
ric AO-type configurations are only stable when C a < 12. 
In the range 12 < C a < 17 the initially purely dipolar 



AO solutions evolve towards a symmetric SO solution af- 
ter about one diffusion time. Furthermore, when C a > 17 
mixed mode solutions are excited. The butterfly diagrams 
for the cases C a < 17 are characterised by strong activity 
at high latitudes, the maximum being near 50°, see Fig.0 
The migration of the activity belts is equatorward at high lat- 
itudes and poleward from the equator up to latitude of about 
45°. This demonstrates that the included surface shear layer, 
leading to equatorward migration in the case of constant in 
radius a-effect, plays a less important role when a concen- 
trates at the bottom of the convection zone. The parameter 
range within which the dipolar AO solution is stable in the 
nonlinear regime is quite narrow. Adding radial and latitudi- 
nal components of the pumping vector or a meridional flow 
pattern is observed to widen this range (see below). 

3.3.1 Effect of the general pumping effect 

The local convection models of KKOS yield all the three 
components of the pumping vector 7. In the present axisym- 
metric investigation, the azimuthal pumping 7^, however, 
would appear only through its radial and latitudinal gradi- 
ents. These terms would act analogously to the differential 
rotation but due to the small value of C a /Ca, these gradi- 
ents are negligible in comparison to the differential rotation. 
Therefore we consider here only the radial and latitudinal 
pumping effects. 

According to the local models, the radial pumping is 
directed downwards throughout the convection zone in the 
slow and moderate rotation regimes, i.e. Co < 4. The radial 
pumping velocity is non-zero also in the case of no rotation 
due to the fact that it is mostly caused by the diamagnetic ef- 
fect (see e.g. Ossendrijver et al. 2002). In the slow and mod- 
erate rotation regime the dependence on rotation (depth) is 
indeed rather weak, although 7,. tends to somewhat increase 
at the poles and decrease at the equator as function of rota- 
tion. The trend is most clear in the rapid rotation regime, 
i.e. Co = 10, where the radial pumping is very small or even 
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positive at low latitudes and at the equator. In the overshoot 
region j r vanishes. We adopt a simplified profile which cap- 
tures the main aspects of the results of KKOS by setting 



7rM) 



1 

'2 a ° 



exp 



tanh ( - 



di 



tanh 



r — ri 



(r - r c ) 2 

4 



sin 9 + 1 



(20) 



where r c = 0.7, n = 0.975, di = 0.015, d 2 = 0.1, and 
c?3 = 0.25. The maximum magnitude of the vertical pump- 
ing is taken to be the same as that of the a-effect, i.e. ao, as 
is suggested by the local results.. 

The latitudinal pumping velocity jg is directed equator- 
wards throughout the convection zone. This effect is rota- 
tionally excited (Krause & Radler 1980), and thus it shows 
a strong rotational dependence so that in the bottom of the 
convection zone it has a peak magnitude of « 2ms" 1 
whilst it is negligible at the surface. Latitudinally ^g is 
strongest near the equator at around 15° with decreasing 
values towards the poles. Similarly to the radial component, 
7e vanishes in the overshoot region. We adopt again a sim- 
plified profile described by 



tanh (V £ )~ tanli ( !1 ^) 



sm 



'cos( 



(21) 

where r 2 — 0.875 and d± = 0.075. The magnitude of the 
pumping effect in the mean-field equations is controlled by 
the nondimensional parameter C 7 ; since the pumping effect 
is basically a part of the a-effect and the maximum magni- 
tudes of both effects were very similar in the study of KKOS 
we set C 1 — C a . 

Adding the pumping velocities defined via Eqs. ( I20t 
and J21t into the model already employing the a-effect 
from convection calculations, widens the region in which 
the dipolar solution is the preferred one (up to C a « 50). 
The dashed line in Fig.[6]indicates that the pumping effects 
make the dynamo easier to excite. This is logical, due to the 
fact that the radial pumping is generally downward, i.e. to- 
wards regions of larger a-effect. The latitudinal pumping in 
the deeper layers has a similar effect, i.e. it transport mag- 
netic fields to lower latitudes. Furthermore, using our stan- 
dard values C a — 15 and Cq = 1.5 ■ 10 4 , the cycle period 
is also significantly longer in this case, i.e. t cyc w 11 years 
instead of approximately seven years in the case without 
the pumping effects. Another clear effect is that the migra- 
tion of the activity belts is now equatorward also at low lat- 
itudes, see Fig. [S] Both of these results seem to be direct 
consequences of the added latitudinal pumping; otherwise 
identical model with jg = yields a butterfly diagram very 
similar to Fig.Q 

3.3.2 Effect of field direction dependent pumping 
effect 

The off-diagonal components of the (symmetric) tensor 
a contribute to the field direction dependent pumping 



(Kitchatinov 1991; Ossendrijver et al. 2002; KKOS). There- 
fore, in the general case, we need to introduce ay compo- 
nents a r g, a r< j), and a^g into the model. However, according 
to the results of KKOS, a^g is clearly the largest in magni- 
tude of the components entering the equation for A. The 
local calculations show that this component peaks strongly 
near the equator at or around latitude 15°, and the maximum 
absolute magnitude is again roughly 2ms" 1 . We introduce 
a profile 



oe<pg{r,9) = -a 



tanh 



di 



tanh 



r — r3 



sin 12 9 cos 2 9 , (22) 

where r% = 0.85 and d 2 = 0.1. The main effect of a^g 
is that it contributes to the radial pumping of the azimuthal 
field via 7^ = j r — a^g, where 7^ is the pumping velocity 
for the azimuthal field. The local results of KKOS indicate 
that the azimuthal field can be pumped upward at latitudes 
|0| < 45°. Due to the narrow latitude and radius range 
where a^g is appreciably large, its effects on the surface 
are hardly visible (compare Figures[8]and[9}. 

3.4 Effect of meridional flow 

Solar surface indicators and local helioseismology (Zhao & 
Kosovichev 2004; Komm et al. 2004) have revealed that 
there is a mean poleward flow of 10-20 ms" 1 near the so- 
lar surface. Although no observational evidence exists as of 
yet, it is generally accepted that there is a return flow at 
some greater depth. However, it is also not known whether 
there is only one cell in radius or more. The simplest case is 
to assume that there is only one cell. We adopt a modified 
version of the flow used by Dudley & James (1989), defined 
via a potential for the poloidal flow 

ip = x 2 sin{-jtx)P%(cos9) , (23) 

where a; = (r — r c )/(R— r c ), and P% (cos#) = i(3cos 2 9 — 
1). Here x describes the effect of decreasing density towards 
the top of the convection zone. The radial and latitudinal 
velocity components are obtained via 



1 



d 



(24) 



ug 



(25) 



89 V 09 J ' 

1 d dip 
r dr 89 

which are written out explicitly in Eqs. iA3\ and JA4I> in 
Appendix A. 

Fixing the maximum value of the flow at the surface to 
10ms" 1 yields C v = 75 or Cu/Cq, = 0.005 for the stan- 
dard model where C a — 15. Introducing this flow into the 
model discussed in Sect. 13.3.21 which already includes both 
pumping effects, we obtain the butterfly diagram shown in 
Fig. The most visible effect of the meridional flow is 
the appearance of the polar branch at latitudes 9 > 60°. 
Furthermore, the equatorward propagating activity belts ap- 
pear now at slightly lower latitudes, although the difference 
to the observed solar activity is still quite large. The criti- 
cal dynamo numbers in the case with meridional flow are in 
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general somewhat higher than in the case with only radial 
and latitudinal pumping taken into account. Furthermore, 
when Cn is smaller than approximately 1.5 • 10 4 , the sym- 
metric SO mode is somewhat easier to excite than the AO 
mode (not shown). 

In the nonlinear regime the solution is dipolar (AO) for 
C a < 15, quadrupolar in the range C a sw 15 . . . 30, and 
mixed parity solutions are obtained when C a > 30 keeping 
the ratio Cy /Cn = 0.005 fixed. If Cu / Cn is incereased the 
activity belts are shifted closer to the equator, but the SO and 
mixed modes become somewhat easier to excite. Increas- 
ing the dynamo numbers is associated with an increasing 
cycle period, but with the present profiles for the a-effect 
and turbulent pumping the solutions in the solar range, i.e. 
t cyc w 22 years, all consist of mixed modes. 

We note here that in the present models the magnetic 
Reynolds number is equal to Cjj, i-e. Rm = Cjj « 0(100). 
This means that we are not in the advection dominated 
regime realised in many flux transport dynamo models (e.g. 
Dikpati & Charbonneau 1999; Bonanno et al. 2002) where 
values in excess of 10 3 are common. The value of Rm, how- 
ever, indicates that the present models are also not fully dif- 
fusion dominated, but rather represent a distributed dynamo. 

3.5 The effect of the surface shear layer 

In order to study the influence of the surface shear layer 
on the solution we take the same setup as in the previous 
section and turn off the shear above r = 0.95i?. The re- 
sulting butterfly diagram is shown in Fig. \^\ Comparing 
this with Fig. 1101 shows some differences, most visibly the 
less extended active regions, and a narrow region of pole- 
ward migration near the equator. In contrast to the case of 
uniformly distributed a-effect (see Sect. 13. 2\ for which case 
the surface shear plays a significant role, the importance of 
the surface shear is very small if the a-effect more concen- 
trated near the bottom of the convection zone. 

3.6 Phase relation of B r B^, 

We note briefly that in our models the correlation of the ra- 
dial and azimuthal magnetic field components, i.e. B r B^, 
is mostly negative (see Fig. ^] for a typical result) in the 
latitude range where the activity belts are located in accor- 
dance with solar observations (e.g. Stix 1976). This relation 
has commonly been used as a constraint for solar dynamo 
models (e.g. Stix 1976; Schlichenmaier & Stix 1995; Bran- 
denburg 2005) although this has recently been criticized by 
Schiissler (2005) who explains this relation as the conse- 
quence of the systematic tilt of the bipolar regions due to 
the action of Coriolis force on rising flux ropes in the con- 
vection zone. 

4 Conclusions 

We have studied kinematic mean-field solar dynamo mod- 
els using helioseismologically determined rotation profile 



and a-effect and turbulent pumping from the local convec- 
tion calculations of Kapyla et al. (2006). We show that if a 
simple a-effect profile proportional to cos 8 in latitude and 
constant in radius is applied, the results are very sensitive to 
the presence of a surface shear layer at r > 0.95i?. Omit- 
ting this feature, the migration of the activity belts is equa- 
torward near the poles, weak at midlatitudes, and poleward 
at low latitudes. The inclusion of the shear layer results in 
equatorward propagation at virtually all latitudes. 

The a-coefficients from local convection calculations 
tend to be much more concentrated near the bottom of the 
convection zone and also closer to the equator with peak 
values occuring at around latitude 30°. In this case, how- 
ever, the critical dynamo numbers are significantly larger 
due to the more concentrated nature of the a-effect. The 
solution is dipolar when C a < 12 and quadrupolar (SO) 
if C a is larger even if a purely dipolar initial condition is 
used. Adding a latitudinal pumping velocity of the order 
of 2ms" 1 at the bottom of the convection zone remedies 
this problem and helps to shift the activity belts closer to 
the equator, whereas the vertical pumping seems to have a 
much smaller overall signifigance for the stability and over- 
all structure of the solution. The latitudinal pumping also 
helps to restore the equatorward migration of activity belts, 
lost by the implementation of the a-effect from convection 
calculations. Adding the pumping effects also makes the cy- 
cle period significantly longer. Additional upward pumping 
of the azimuthal field component through the component 
a^e at near equator regions has very little effect on the so- 
lution. In all the models investigated, the correlation BrB^ 
is mostly negative, which is in agreement with the observa- 
tions, but none of the models produces poleward migration 
at high latitudes (polar branch). 

If a single cell meridional flow patterns where the flow 
is of the order of 10 ms -1 and poleward near the surface 
and significantly smaller and equatorward at larger depths 
helps to shift the activity belts further closer to the equator, 
although they still appear too high (|0| w 5° . . . 60°) in 
comparison to the Sun. A region of poleward migration also 
appears at high latitudes. If the upper shear layer is removed 
from this model, the butterfly diagram changes only little in 
contrast to the models where there was an appreciably large 
a-effect also near the surface. 

The distributed dynamo model investigated in this paper 
correctly reproduces many of the general features of the so- 
lar cyclic activity, including realistic migration patterns and 
correct phase relation, but some problems persist. These in- 
clude the excitation of the dipolar (AO) solutions only within 
a narrow range of parameters with the anisotropic a-effect 
combined with meridional circulation, activity at too high 
latitudes for all models investigated, and the somewhat too 
short dynamo cycles. With further tuning of the various pa- 
rameters going into the model (magnitudes of a- and j- 
effects, strength and profile of the meridional circulation) 
it would most likely be possible to obtain a solution repro- 
ducing the missing features, as well. Instead of going fur- 
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ther into the domain of tuning in the present study, we note 
that further work in the context of local convection mod- 
elling is needed in order to obtain the totally missing infor- 
mation on the anisotropics of turbulent diffusivity (work in 
progress). It is also necessary to consider the nonlinear sat- 
uration process of the a-effect in more detail, for instance 
in the form of a dynamical a-effect with helicity fluxes (e.g. 
Brandenburg & Subramanian 2005). Further improvements 
of the present mean-field model include the relaxation of 
the kinematic approach to include full dynamics, for which 
purpose information from local convection calculations also 
exist (e.g. Kapyla, Korpi & Tuominen 2004). With the inclu- 
sion of hydro- and thermodynamics the nonlinear feedbacks 
of the kind investigated by e.g. Rempel (2006) could be in- 
cluded in a self-consistent way. 

Despite of the shortcomings discussed above, we feel 
that the simple model presented in this paper demonstrates 
two important aspects concerning solar dynamo theory. 
Firstly, as recently discussed at length by Brandenburg 
(2005) the surface shear layer may be important in shap- 
ing the solar dynamo; according to our results this is indeed 
to be expected if the a-effect has an appreciable magnitude 
near the surface. Secondly, we are able to demonstrate that a 
distributed ail-dynamo model applying a realistic rotation 
profile and meridional flow can reproduce the correct phase 
relation and migration of the activity belts. Although merid- 
ional circulation is an essential ingredient also in our model, 
the magnitude of this effect is of the same order as the tur- 
bulent inductive effects and the turbulent diffusivity is an 
order of magnitude larger than in the advection dominated 
flux-transport dynamos. 
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A Equations 

The equation for A reads 



9 A .„ „ TT J A dA 



, rr . ( cot 9 . 1 dA 

+ C Q a 0r ~CuU e )[ — 

\ r r d6 

+ CaCtttxftB 

fp 2 A 2 OA 1 d 2 A | cotOdA 

T ^ t I dr 2 r dr r 2 09 2 r 2 89 



(Al) 



where Ui = Ur r + 7&7i- 



dB 

a7 



Cnr sin 9 



Cn sin 9 



Cu 

0Q /cot0 



Or ( 

on [A 0A\ 
~09 \~r + ~0r) 



1 QA\ 



A + , , 
r r 09 I 



Cu M- + C u -^(cot9A+ , 
r r Or \ 09 J 



) 



Or 



\ r 



0A\ 



© 0000 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhdm 



www.an-journal.org 



Astron. Nachr. / AN (0000) 



9 




Fig. 1 Upper panel: rotation profile adapted from the he- 
lioseismic results. The lower panel shows the angular ve- 
locity O normalised by the core angular velocity, S! C orc = 
430 nHz, at latitudes 0° (equator, thick solid line), 15° 
(dashed), 30° (dot-dashed), 45° (triple-dot-dashed), 60° 
(solid), 75° (long-dashed), and 90° (pole, dotted line). The 
dashed vertical lines at r — 0.70i? and r = 0.95i? de- 
note the bottom of the convection zone and the depth above 
which the surface shear layer is located, respectively. The 
thin horizontal lines above r = 0.95R give £1 in the case 
where the surface shear is neglected. 
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Fig. 2 Critical dynamo numbers C" 1 ' as function of Cn 
for the a-profile given by Eq. dl6l . The solid and dashed 
lines give the results for the rotation profile including the 
surface shear, shown in Fig. Q an d using the aQ and 
a 2 Vi approximations, respectively. The dot-dashed (afl) 
and triple-dot-dashed (a 2 £l) curves give the corresponding 
results with a rotation profile where the surface shear layer 
has been removed. The thick solid line describes the Sun, 
i.e. where C q /Cq = 10~ 3 . 
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Fig. 3 Butterfly diagram of the toroidal magnetic field 
from r = 0.99R from a model with C a = 15 and Cq — 
1.5 • 10 4 with the rotation profile shown in Fig. Q with the 
surface shear at r > 0.95R artificially turned off. 
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Fig. 6 Critical dynamo number C£ rlt as a function of Cn 
for the antisymmetric AO mode with the a-profile from 
convection calculations, see Fig. [5] (solid line). The addi- 
tional curves are with latitudinal and radial pumping ac- 
cording to Eqs. (I20> and \2\\ (dashed line), and with merid- 
ional flow according to Eqs. (IA3I and \AA\ (dot-dashed 
line). The thick solid line describes the Sun, i.e. where 
C a /C n = 10- 3 . 
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Fig. 7 Butterfly diagram of the toroidal magnetic field 
from r = 0.99R from a model with C a = 15 and Cq 
1 .5- 1 '. The cycle period is approximately p Vi . w 7.0 years, 
a^^ -profile as shown in Fig. |5] 



90 



4 



-4 




Fig. 10 Same as Fig. [9] but with added meridional flow 
with max(iie) = 10ms _1 at the surface, with which choice 
Cjj = 75. The cycle period is about 1 1.5 years. 
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Fig. 8 Same as Fig. [7]but with added radial and latitudinal Fig. 11 Same as Fig. 1101 where the surface shear above 
pumping velocities as defined via Eqs. J20l > and i2l\ . The r — 0.95R has been artificially turned off. The cycle period 



cycle period is t cyc w 12 years. 



is about 9.5 years. 
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Fig. 9 Same as Fig.|8]but with added field direction de- 
pendent pumping according to Eq. J22b . The cycle period is 
about 12 years. 
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Fig. 12 The correlation B r B^, from the model shown in 
Fig- E| Black denotes negative and white positive values. 
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